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1.  INTRODUCTION 


We  report  here  the  results  of  a  study  of  the  feasibility  of  measuring 
crosswinds  by  observing  fluctuations  in  the  image  formed  by  a  lens  observing 
an  arbitrary  scene.  The  study  is  especially  concerned  with  observations  in 
the  10  ym  wavelength  region  used  by  the  FLIR  thermal  sight  system,  and  how 
this  passive  system  might  be  used  to  obtain  crosswind  information.  This 
report  is  submitted  to  fulfill  the  requirements  of  Project  Order  ASL798003 
issued  by  the  U.  S.  Army  Electronics  Command,  Atmospheric  Sciences  Labora¬ 
tory,  White  Sands  Missile  Range,  New  Mexico. 

A  previous  report*  discussed  the  feasibility  of  measuring  crosswinds 
by  observing  the  scintillation  of  the  natural  background  radiation.  The 
analysis  technique  was  that  used  for  previous  passive  wind  sensing  devices 
operating  at  visible  wavelengths.  A  new  theoretical  development  is  required 
to  understand  the  origin  of  image  plane  fluctuations.  The  theory,  presented 
in  the  next  section,  is  compared  quantitatively  with  data  obtained  from  two 
experiments.  Additional  quantitative  comparisons  which  involve  extensive 
computer  calculations  will  be  available  in  a  future  publication. 

In  one  experiment,  we  made  continuous  observations  of  crosswinds  at 

visible  wavelengths  in  real  time;  in  the  other,  tape  recorded  electrical 
signals  from  the  FLIR  instrument  were  analyzed.  Both  experiments  were 
conducted  over  the  same  500  m  test  range.  The  visible  wavelength  experiment 
was  conducted  under  more  favorable  conditions  than  the  FLIR  experiment  and  the 
simpler  test  arrangement  was  easier  to  analyze.  The  FLIR  analysis  was 
complicated  because  it  uses  a  system  with  intermittent  sampling,  and  has 
other  problems,  but  actual  use  of  the  instrument  was  extremely  helpful  in 
understanding  the  practical  difficulties  of  adding  an  optical  crosswind 
capability  to  this  instrument. 


2.  THEORY  OF  WIND  MEASUREMENT  BY  OBSERVATION  OF  IMAGE  FLUCTUATIONS 


We  can  analyze  the  FLIR  problem  using  the  classical  Huygen- Fresnel 
(H-F)  approach:  Figure  1  illustrates  a  scene  u^pp  a  distance  away  from  a 
lens  which  images  the  scene  a  distance  f  +  5  in  back  of  the  lens.  The  image 


u.(p)  is  given  by 
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The  field  u„(P-,)  is  related  to  the  field  incident  upon  the  lens  u7(p,,) 
by  the  lens  transfer  function 


u3(p2)  =  u2(p2)  e‘lk°2  /2f  W(p2)  •  (2) 

where  W  is  the  aperture  weighting  function,  assumed  Gaussian,  W(p2)  = 
exp  [r‘‘/2rQ2]  where  rQ  is  (approximately)  the  lens  radius,  lo  complete  the 
equation  for  the  imaging  of  the  actual  scene,  wc  need  the  expression  for 

O 

u,(p7)  in  terms  of  the  scene  itself.  Applying  the  extended  (H-F)  approach", 
•o  nave 
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-ik  +  ikz 


2 

,  -i  k  (P7  .  Pi)  +  1  kPi  /(2z)  U»fo  o-o  z  tl 

/  ,rP  el-2-1  1  u(Pj) 


(3) 


where  z  is  the  longitudinal  coordinate. 

Equation  (3)  contains  most  of  the  important  physics  of  the  problem.  This 
equation  differs  from  the  (HF)  integral  Eq.  (1)  by  the  factor  exp  M  which 
contains  the  effects  of  the  refractive  turbulence  intervening  between 
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source  and  receiver.  The  HF  approach  assumes  that  you  can  break  up  any  field 
distribution  e.g.,  U1(P1)  in  Eq. (3)  into  wavelets  originating  at  each  point 
pj,  propagate  them  to  the  receiver  and  recombine  the  effects  of  each  wavelet 
at  the  receiving  plane  to  obtain  the  total  field.  The  addition  of  the  factor 
exp  [i'j]  accounts  for  the  effects  of  refractive  turbulence  on  each  wavelet  as 
it  propagates  from  its  origin  at  (p^,0)  to  the  receiver  at  (p2>z).  The 
function  \jj  =  y  +  is.,  wheTe  y  is  the  logamplitude  and  s1  is  the  phase,  is  the 
Rytov  solution  '  for  the  logarithmic  field  of  a  spherical  wave  propagating 
through  refractive  turbulence. 


Inserting  Eqs.  (21  and  (3)  into  Eq.  (1),  we  have  the  expression  for  the 
image  field  in  the  presence  of  turbulence  in  the  form. 
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The  image  plane  of  a  lens  is  defined  by  the  relation,  z  1  +  (f  +  £}  ^  - 
=  0  or  where  6  ~  f  (f/zl .  Letting  d  =  f  +  6  in  the  image  plane  we  have 
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Note  that  in  the  absence  of  refractive  turbulence  (ij>  -*■  0)  and  an  infinite  lens 
(W  -*■  1) ,  we  obtain 


u. (p)  «  M'1u(-M-1p),  (6) 

X  ~  <v 

the  reflected,  magnified  (M  =  d/zl  image  that  we  expect  from  geometric  optics. 
The  effect  of  refractive  turbulence  will  be  to  defocus  the  image  and  cause  it 
to  dance  around  rapidly  in  the  image  plane. 


To  measure  transverse  horizontal  wind  speed,  we  detect  the  scintilla¬ 
tions  with  two  detectors  of  radius  aQ  separated  in  the  horizontal  direction 
an  amount  l/2p^  on  each  side  of  the  beam  axis.  The  detected  signal  takes  the 
form 
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* 

where  I(p)  =  u. (p)  u.  (p) .  By  calculating  the  time-lagged  covariance  of  the 
~  2 

two  signals  Cj/<I0>  given  by 


CT(pt,T)  =  <I(pt/2,t)  I(-pt/2,  t+T)> 


<Xo> 


(8) 


where  the  angle  brackets  indicate  an  ensemble  average,  <IQ>  is  the  mean 

irradiance  in  the  weak  turbulence  limit  and  statistical  stationarity  is 

assumed.  We  can  then  compute  the  slope  of  this  quantity  m  =  CT(p. ,t) 

o  T  I  —t. 

which  will  be  proportional  to  a  weighted,  path-averaged  transverse  wind 
velocity^.  Substituting  Eq.  (5)  into  Eq.  (7)  and  performing  the  p 
integration  using  standard  tables,  we  obtain 
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where  Kfp^)  =  exp  f - /2t ^  ].  Ohviously  it  is  extremely  cumbersome  to  write 
Cj(pt,Tj  directly  so  we  will  outline  some  simplifying  arguments  first. 


Cj  contains  two  factors  of  interest: 
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(Subtracting  unity  in  Eq.  (11)  is  valid  only  in  the  weak  turbulence  limit,  which 
we  assume  in  the  calculation  below.)  The  quantity  II  is  the  fourth  moment  of 
the  spatial  field  distribution.  If  the  u's  are  representable  as  joint  Gaussian 
spatially  ^coherent  random  fields,  then^ 
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This  assumption  simplifies  C^.  to  the  form 
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the  next  problem  is  to  evaluate  the  P's.  If  we  make  the  reasonable  assumption  that 
the  bulk  of  the  propagation  effects  arise  from  the  influence  of  phase  fluctua¬ 
tions,  Lee  et  al.2  have  shown  that 
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P  +  1  =  exp{-  ^  [DfPj-Pj  ,  £2"£2  ,0^  '  D^£l'£l  '  £2"£2  ,T-* 
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where  D  Is  the  two  point  source,  two  receiver  "wave"  (or  approximately  "phase") 
structure  function  defined  by  D(x,y,T)  -  <  (x  *  ,y '  ,x ' )  +  ^  (x’+x,  y'+y,  t'+t)]^ 
and  equal  to 


nrx,y,T)  =  4-rrk2z  /  dt  /  d2K  4>n(K,t) 


where  is  the  refractive- index  power  spectral  density,  t  is  the  normalized 
path  position,  z  is  the  total  path  length  and  v  is  the  wind  velocity.  To 
further  simplify  the  calculation,  we  make  the  weak  turbulence  assumption, 
valid  far  TR,  exp[A]  1  +  A.  Further,  we  drop  terms  that  are  time  indepen¬ 
dent  and  therefore  do  not  affect  the  slope  at  zero  time  lap.  Under  these 
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approximations,  P  and  P~  could  be  obtained  by  letting  p  =  p  ,  p  - 
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To  simplify  Eq.  (13),  we  insert  Eq.  (13)  and  Eq.  (16)  and  make  the 
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changes  of  variables  £  =  ,  2n  =  =  P2  -  P2  ’  2Dl  =  ~2  +  ~2  ’ 
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C2  =  P2  -  p  ,  2n2  =  P2  +  P2  and  integrate  over  rj1  and  n2  to  obtain 
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valid  for  ka  /d  >>  r  .  (For  our  problem  a  /d  ^  10  k  ^  2tt  x  10^  m  and 
00  0 

r  =  7  x  10-2  m  so  ka  /d  ^  630  m"1  and  r  _1  ^  14  m-1.)  To  continue,  we  let 
o  o  o 

2y  =  ^  ,  x  =  -  ?2  and  perform  the  x,  y  integrations,  resulting  in 

the  expression 
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where  p  *  p  z/d,  a  =  a  z/d  and 
~s  ~t  so 
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Now  we  must  insert  some  model  of  our  FL1R  target.  A  useful  model  and  one 
directly  applicable  to  our  experiments  is  the  horizon  where 
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Note  for  this  form  of  I  that 
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Inserting  this  into  Eq.  (18)  we  can  start  the  tedious  but  straightforward 

process  of  integrating  over  5  and  n'.  Performing  the  calculation, 

for  weak  refractive  turbulence  and  for  z/(k  r  tx  )  «  1,  we  obtain 
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where  we  have  normalized  C_  by  <1  >  .  (Note  the  last  two  terms  in  Eq.  (18) 

X  o 

contribute  nothing  to  the  slope  at  zero  lag  for  the  particular  target  that 
we  have  assumed.) 


From  Eq.  (21)  we  compute  the  slope  at  zero  lag  m  contributed  by  the 
x-directed  or  horizontal  wind  component  (x  ||  p^). 


m  =  /  dt  C  ^(t)  v  (t)  W(t)  , 
n  x 

o 


(22) 


where  the  wind-C  weighting  function  W  is  given  by 
n 


W(t)  =  0. 528irk2z  /  dK^K^  sin  {Kj>o(l-t)]  /  dK  K~'1/3 

v  v  _ oc  y 


X  X 


X  s 


2  2  2  2  2  2 

-K  t  r  12  -  K  (1-t)  a  /2 

o  s  .  2 

e  sin 


fv2  /.  v  r  o  /iK  (l-t)u  \1 

— s)J 


(2  i) 


(Note  the- y-directed,  vertical  wind  component  will  not  contribute  to  m.)  We  have 

3  4 

inserted  the  Kolmogorov  form  ’  of  the  refractive  index  spectrum  in  Eq .  (23),  that 
is. 


?  -11/3 

$  (K, t)  =  0.033  C  Z(t)  K  (24) 

n  ~  n 

valid  if  L  3  <<  K  &  3  where  L  and  2.  are  respectively,  the  outer  scale 

o  o  o  o 

and  microscale  of  the  refractive  index  fluctuations. 

Although  Eqs.  (21)  and  (23)  were  derived  for  phase-variation-induccd 
image  fluctuations  and  their  ultimate  contribution  to  irradiancc  fluctua¬ 
tions,  the  form  of  these  equations  suggest  a  convenient  analog.  When 

9  . 

comparing  Eq.  (21)  with  Eq.  (9)  of  Clifford  et  al.,  it  is  apparent 
that  a  completely  analogous  calculation  to  the  FLIR  problem  with 
a  horizon  target  is  wind  sensing  using  two,  half-moon,  spatially  incoherent 
sources  with  their  flat  edges  oriented  in  the  x  direction  and  displaced  an 
amount  pg  parallel  to  their  edges.  In  this  analog,  these  sources,  that  are 
located  at  the  FLIR  target  plane,  must  be  coded  and  observed  with  two 
co-located  receivers  at  the  aperture  of  the  FLIR.  The  equivalent  sources 
and  receivers  must  have  the  amplitude  weighting  function 


15 


-i  2  2  2  2  2 

exp[-(x~+y  )/as  ]  for  y  £  0  and  exp[-(x  +y  )/r  ]  , 

9 

respectively.  We  have  previously  analyzed  this  problem  and  the  applicable 
wind  weighting  functions  peaked  closer  to  the  receiver  (FLIR) ,  than  toward 
the  transmitter  (target)  by  an  amount  depending  on  the  ratio  of  the  source 
.u.d  ceiver  radii.  As  will  be  shown  in  the  next  section,  this  is  quali¬ 
tatively  observed  in  our  FLIR  experiments.  A  computer  analysis  of  Lq.  (23) 
will  be  undertaken  to  give  precise  results  for  our  experimental  conditions. 

A  rough  idea  of  the  shape  of  the  wind  weighting  function  and  its  de¬ 
pendence  on  (r  /a  )  can  be  seen  as  follows:  If  we  note  for  r  <  a  that  Kt  < 
1  .  o  s  7  o  s  — 

r  in  the  integrand  of  Eq.  (23)  then,  in  the  sin'G  term,  if  we  ignore  the 

°  i  2  2  2 

end  points  where  K(t)  =  0  and  use  sin  t  ^  2,  e  =  z/(krQ  )  <  1  and  sin  6  a-  0 

is  a  reasonable  approximation.  With  this  assumption  and  a  change  of  vari¬ 
ables  to  K'  =  K(l-t),  we  obtain 


W(t)a  - jtj- 

a-ty'* 


/  dKY\'  sin(K/P<)  /  dKv'(Kv'2+K^) 


2,1/6 


,2‘2-2-*  -K2a2/2 


exp 


r  '  £.1  £n 
-K  t  r 
_ o 

- 2(l-t) 2  . 


e  S  '  [1  -  erf2(iKy’as/2)  . 


r:s) 


We  note  that,  if  p  a  ,  K'  /2/a  in  the  important  region  of  integration 

i  ?  *  s 

in  noth  K  and  K  ,  then  we  obtain  as  a  rough  approximation 
y  ’ 


W (t)  a  - Wrr  e 


(1-t) 


4/3 


(26) 


A  plot  of  this  function  is  shown  in  Fig.  2  for  r0/<*s  =  1  (dotted  line)  and 
ro/ag  =0.5  (solid  line).  The  plot  compares  favorably  with  the  experimental 
results  in  Section  3.  In  general,  the  weighting  function  Eq.  (26)  and  the 
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approximations  above  indicate  that  the  peak  of  the  weighting  function  is  at 

t  ■  (1  +  r  /a  and  the  most  important  scale  size  for  producing  fluctu- 
m  os 

ations  is  ?.  ^  /2u  (r  *  +  a 

o  s 

3.  EXPERIMENTAL  MEASUREMENTS 

3.1  Wind  Measurements  at  Visible  Wavelengths 

All  of  the  field  measurements  for  this  study  were  made  at  Table  Mountain, 

a  mesa  having  a  flat  area  of  about  1.5  by  3  km,  and  located  12  km  north  of 

Boulder,  Colorado.  Figure  3  shows  a  view  lookirg  down  the  500  m  test  range 

used  for  the  experiments.  The  two  trailers  and  the  surface  of  the  mesa 

surrounding  them  comprised  the  scenes  observed. 

For  the  visible  experiments,  a  Schmidt-Cassegrain  telescope  of  22.3  cm 
diameter  was  focused  to  form  an  image  on  a  dual  photodiode  by  eyepiece 
projection,  as  in  Fig.  4.  The  angular  fields  of  view  were  determined  by  the 
size  of  the  image  on  the  photodiode  surfaces.  The  electrical  signals  from 
the  photodiodes  were  then  processed  as  follows:  The  normalized  covariance 
versus  time  delay  was  observed  in  real  time  by  a  correlator.  Two  particular 
•iroporties  of  the  covariance,  its  slope  at  zero  time  delay  and  an  analysis 
ha  ;od  upon  a  sampling  of  the  function  at  14  time  delays'^  were  computed. 

Av>  -aged  values  of  these  quantities,  together  with  information  from  a  pro- 

2  .  • 

peller  anemometer  array  and  an  optical  C  meter,  were  logged  bv  a  mini- 

n 

computer  which  also  computed  the  weighting  functions  for  the  optical  wind 
measurement . 

Typical  covariance  functions  are  shown  in  Fig.  5.  According  to  the 
theory,  the  propagation  characteristics  for  the  optical  geometry  used 
should  >e  very  similar  to  that  which  would  exist  if  the  system  were  replaced 
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jaindwoo 


Figure  4.  Block  diagram  of  the  experiment  at  visual  wavelengths. 


by  a  detecting  aperture  the  size  of  the  receiver  (22.3  cm)  and  two  radiating, 
halt'-moon,  apertures  the  diameter  of  the  projected  area  of  the  detectors 
on  the  scene  (40  cm) .  With  this  geometry  one  would  expect  from  the  theory 
a  peak  response  to  winds  about  t  a,  1.22  rQ/ois  =  0.68  of  the  way  from  the 
scene,  and  maximum  sensitivity  to  irregularities  of  spatial  wavelength  about 
(2ttK~  1 )  %  2-rras/(3/2)  =  30  cm.  Since  the  average  crosswind  was  measured  by 
anemometers  during  the  covariance  function  measurement,  the  spatial  wave¬ 
length  can  be  estimated.  It  turns  out  to  be  30  cm,  in  agreement  with  the 
prediction.  The  weighting  function  was  also  measured  (Fig.  6)  and  its 
peak  response  was  as  predicted. 

In  the  experiment,  we  looked  at  the  horizontal  edge  of  the  window  of 
the  trailer  which  was  dark  against  the  light  outside  of  the  trailer.  The 
horizontal  edge  also  assured  that  both  detectors  saw  the  same  scene  spatial 
spectrum.  The  theory  of  course  shows  that  any  measure  of  wind  by  this 
technique  is  a  function  of  the  scene  spatial  spectrum,  and  since  the 
detectors  are  observing  relatively  small  portions  of  the  scene,  in  general 
wide  variations  in  the  spectrum  may  be  expected.  In  practice  this  means 
that  wind  information  can  only  be  obtained  from  certain  portions  of  the 
scene.  At  first  it  was  thought  that  by  averaging  the  information  from 
many  pairs  of  detectors  observing  the  scene,  the  signal -to-noi se  of  the 
wind  measurement  could  be  improved  by  something  like  the  square  root  of 
the  number  of  pairs  of  detectors.  Experience  has  shown  that  this  would 
not  be  the  case.  In  fact,  one  pair  of  detectors  at  a  particularly  favorable 
location  in  the  scene  will  very  likely  have  better  signal-to-noise  than 
the  average  of  many  pairs.  With  many  pairs  of  detectors  available, 
h  v,  aver,  one  would  have  the  advantage  of  being  able  to  choose,  by  criteria 
such  as  cintillation  level  and  correlation,  the  detector  pairs  having 
the  best  signal-to-noise  for  wind  measurement. 

Figure  7  shows  comparisons  of  wind  measurements  obtained  by  the  slope 
and  covariance  techniques  with  the  average  of  10  propeller  anemometers  along 

the  light  path  for  a  sunny  period  of  3  hours  in  the  afternoon  of  October  3, 

2 

1979.  An  optical  measurement  of  Cn  is  also  shown.  The  optical  measurements 
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Figure  7.  Optically  derived  wind  speeds  c 
10  propeller  anemometers  along  the  op 
level  during  the  experiment  is  also  s 


24 


correlate  well  with  the  anemometer  averages;  however,  the  conditions  of  high 
2 

C  and  a  carefully  selected  target  were  ideal  and  it  is  clear  that  the 
system  could  not  be  depended  on  to  operate  in  more  adverse  conditions. 


3.2  Wind  Measurements  from  FLIR  Thermal  Sight  Signals 

Personnel  from  the  Atmospheric  Sciences  Laboratory,  White  Sands  Missile 
Range,  New  Mexico,  provided  an  AN/VSG  Thermal  Sight  of  the  tank-mounted  type, 
set  up  the  system  on  Table  Mountain  and  explained  its  operation.  They  also 
provided  a  high-speed  14- channel  instrumentation  tape  recorder  and  worked 
with  NOAA  personnel  to  record  the  data  which  is  analyzed  and  presented  in 
this  section. 

References  11  and  12  describe  the  thermal  sight.  It  is  described  as  an 
integrated  dav/night  fire  control  periscope  for  the  M60A1(PI)  main  battle 
tank.  It  combines  a  daylight  optical  viewing  system,  a  self-contained  ballistic 
reticle  projector,  and  a  high-resolution  passive  thermal  imaging  system  that 
provides  both  surveillance  and  fire  control  capability  to  commander  and 
gunner  during  day  and  night  battlefield  conditions.  Some  of  the  data  per¬ 
tinent  to  our  study  are  listed  in  table  1.  The  image  is  formed  by  scanning 
a  vertical  array  of  120  IR  detectors  horizontally  across  the  image  plane  with 
a  scanning  mirror,  then  shifting  vertically  by  one-half  the  detector  spacing  for  an 
interlaced  return  trace,  at  the  rate  of  42  of  these  cycles  per  second.  The 
electrical  signals  from  these  detectors  drive  LEDs  which  arc  scanned  with  the 
same  mechanical  assembly  to  form  a  visual  image.  The  video  signals  from  all 
120  detectors,  as  well  as  a  sync  pulse  at  84  Hz,  are  readily  available.  The 
test  setup  for  recording  video  data  from  the  thermal  sight  is  shown  in  Fig. 

8.  Three  adjacent  detector  scans,  the  sync  pulse,  voice  annotation,  and  9 

propeller  anemometers  were  tape  recorded.  A  separate  chart  recorder  con- 

2 

tinuously  recorded  levels  from  an  optical  meter.  The  video  bandwidth  of  the 
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TABLE  1.  THERMAL  CHANNEL  PARAMETERS  FOR  NARROW  FIELD  OF  VIEW 


Field  of  view  (degrees) 
Frame  rate 
Interlace  ratio 
Scan  efficiency 
Germanium  primary  lens 
Diameter 

Effective  focal  length 
Detector 
Type 

Number  of  elements 
Spectral  wavelength 


2.58  x  5 
42  Hz 
2:1 
0.70 

15.2  cm 
27.1  cm 

Mercury  cadmium  telluride 
120 

7.6  to  11.75  micrometers 


FLIR  is  70  kHz,  but  it  was  reduced  to  23  kHz  in  the  recording  process. 

This  bandwidth  resulted  in  a  horizontal  angular  resolution  of  about  60  yrad. 

As  in  the  experiment  at  visible  wavelengths,  we  wish  to  examine  the 
changes  in  brightness  at  a  particular  location  in  the  image  plane,  and 
correlate  these  changes  with  those  at  another  location,  horizontally 
displayed  from  the  first.  Ideally,  one  would  like  to  do  this  for 
all  possible  pairs  of  horizontally  displaced  points.  Two  approaches  were 
considered  for  this  analysis.  One  was  to  digitize  a  small  portion  of 
the  video  record  at  a  very  high  rate  (>50  kHz).  Due  to  the  low  frequencies 
present  in  the  signal,  however,  this  sample  would  have  had  to  be  about 
5  to  10  seconds  long.  Another  way  would  be  to  synchronize  the  beginning 
of  the  digitizing  with  the  video  sync  pulse,  and  then  digitize  a  portion 
of  every  other  trace,  starting  at  the  same  location  on  each  trace.  In 
fact,  portions  of  the  video  tape  recordings  were  digitized  in  the  latter 


way,  at  an  8  kHz  rate.  It  soon  became  clear,  however,  that  this  rate 
would  not  really  be  adequate,  and  that  in  addition,  it  would  not  be  as 
advantageous  to  average  the  correlation  functions  of  many  pairs  of  points 
as  we  first  thought  it  would  be. 

We  had  planned  to  use  another  method  to  preview  the  tape  recorded 
information,  but  it  turned  out  to  be  the  best  method  for  the  entire  analysis. 

Only  one  pair  of  points  could  be  analyzed  at  once,  but  there  was  no  bandwidth 
limitation,  the  traces  could  be  monitored  in  real  time  on  an  oscilloscope, 
there  was  a  high  degree  of  flexibility,  and  most  important,  nearly  instan¬ 
taneous  turnaround  time.  A  block  diagram  of  the  analysis  technique  is  shown 
in  Fig.  9.  The  sync  signal  and  the  video  signal  from  one  of  the  detectors  enters 
the  sample  and  hold  circuit  enclosed  by  the  dotted  line.  This  specially 
designed  circuit  operates  as  shown  in  Fig.  10.  A  time  delay  i  from  the  end 
of  the  sync  pulse  is  adjusted  to  sample  at  the  desired  location  on  the  trace. 

At  this  time  the  signal  amplitude  is  sampled  and  held  until  the  interlace 
sweep  is  completed  and  the  original  trace  is  beginning  again.  The  signal  is 
again  sampled  at  the  same  location  on  the  trace  and  held.  Disregarding 
jitter  and  other  system  noise,  the  fluctuations  in  the  level  of  this  sampled 
waveform  represent  changes  in  brightness  at  a  fixed  location  in  the  image 
plane.  A  second  sampled  waveform  is  generated  by  sampling  and  holding  the 
signal  amplitude  of  the  video  trace  at  a  time  delayed  by  A£  from  the  sample 
times  of  the  first  trace.  The  normalized  covariance  of  these  signals  can 
then  be  computed  in  the  same  way  that  it  was  done  at  visible  wavelengths  for 
the  continuous  signals  obtained  from  fixed  detectors,  horizontally  separated 
in  the  focal  plane. 

Since  the  sampling  rate  is  42  per  second,  signal  fluctuations  above  21 
Hz  cannot  be  observed.  In  fact  because  of  noise  the  frequency  limitation  is 
even  lower;  this  limits  the  maximum  wind  speed  that  can  be  measured,  and  the 
data  shows  this  effect. 
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R  analysis  techniq 


About  75  minutes  of  FLIR  data  were  recorded  on  10  rolls  of  c 

tape  on  February  7,  8,  9,  12,  13,  14,  15,  and  21,  1979.  We  examined  all 
of  the  tapes  but  chose  to  analyze  intensively  tapes  7,  8,  9,  and  10  taken 
on  February  14  and  21.  Covariance  functions  for  2  sample  points  taken 
100  ys  apart  from  data  on  tapes  7,  8  and  9  are  shown  in  Figs.  11  through  13. 
Full  scale  plus  and  minus  on  the  ordinate  represents  a  normalized  covariance 
of  +1  and  -1  respectively.  On  the  abscissa,  the  time  scale  is  30  ms  per  bit 
or  .96  s  full  scale  for  all  figures  except  Fig.  16,  which  is  60  ms  per 
bit  or  1.92  s  full  scale.  The  time  delay  A£  between  the  samples  from 
which  the  covariance  is  computed  is  100  ys,  except  for  Figs.  14,  15,  and 
16,  where  it  is  25,  50,  and  200  ys,  respectively.  The  covariance  estimate 
is  made  from  20  s  of  data  in  each  case.  The  time  sequence  for  the  co- 
variance  functions  progresses  from  upper  left  to  lower  right. 

Generally,  a  covariance  function  with  a  time-delayed  peak  corresponding 
to  the  existing  crosswind  was  observed.  The  series  in  Fig.  12  is  rep¬ 
resentative.  In  this  case,  the  direction  of  the  crosswind  changed  during 
the  run,  shifting  the  peak  from  positive  to  negative  delays. 

Knowing  the  average  wind  speed  from  anemometer  readings  (0.6  m/s)  and 
obtaining  from  Fig.  12  the  average  time  delay  to  the  peak  of  the  covariance 
function  (0.2s),  we  can  estimate  the  separation  in  space  of  the  light  beams 
as  0.6  m/s  x  0.2  s  =  0.12  m.  If  the  sampled  points  on  the  trace  are  sepa¬ 
rated  by  100  ps,  the  angular  divergence  a  of  the  rays  from  these  points  to  the 
scene  is 


a  = 


iuu  us 


10  ms  x  57. 3 9 /rad 


=  8. 15  X  IU 


Roughly  then,  the  wind  weighting  function  should  peak  up  at  a  distance 
(1-t)  from  the  receiver,  or 

0  12 

(1-t)  =  - — - j  =  137  m, 

8.72  x  10'* 

which  is  very  close  to  the  prediction  of  Eq.  (26),  using  a  15  cm  receiving 
aperture  and  detectors  whose  effective  projection  size  at  500  m  are  0.5  m. 
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Using  this  calibration,  we  plot  in  Fig.  18  a  comparison  of  the  propeller 
anemometer  average  versus  the  optical  measurement  of  wind.  While  the 
data  from  tape  8  compares  quite  well,  the  optical  measurement  is  too  low 
for  the  data  on  tapes  7  and  9.  This  result  is  expected  since  the  42  Hz 
sampling  rate  would  limit  the  frequency  response  so  that  for  this  geometry, 

2  m/s  winds  would  be  about  the  maximum  speed  that  could  be  measured. 

Higher  wind  speeds  could  in  principle  be  measured  by  shortening  the 
sample  interval  time  AJl.  Both  the  resolution  of  the  system  and  the  atmosphere 
limit  the  minimum  time,  however.  Figures  14,  15,  11,  and  16  show  the  same 
data  on  tape  8  analyzed  using  sampling  delays  of  25,  50,  100  and  200  ys. 

The  25  ys  data  show  some  time  delayed  covariance,  but  mostly  peak  up 
around  zero  time  delay.  This  is  a  result  of  insufficient  video  bandwidth. 
While  the  FLIR  system  had  sufficient  bandwidth,  the  recording  system  had 
less  than  half  of  the  frequency  response  required  to  record  this  signal 
accurately.  There  is  nearly  enough  bandwidth  available  for  the  measurements 
made  at  50  ys  intervals  (Fig.  15) .  One  can  see  that  the  peak  delay  is 
about  half  of  that  at  100  ys  (Fig.  11)  although  there  is  some  tendency 
to  peak  around  zero  delay.  From  the  theory  one  would  expect  that  some 
peaking  symmetric  with  zero  time  delay  could  occur  due  to  back  and  forth 
tilts  of  the  optical  wave  front  near  the  receiving  aperture.  At  200  ys 
(Fig.  16,  note  time  scale  change),  the  peak  delays  are  twice  as  long, 
as  one  would  expect,  but  also  there  is  evidence  of  less  covariance  at  the 
wider  spacing. 

An  example  of  the  sensitivity  of  the  system  to  a  particular  portion 
of  the  scene  is  shown  by  the  bottom  graph  in  Fig.  12.  It  shows  the  co- 
variance  function  obtained  from  signals  obtained  during  the  return  trace. 

The  field  of  view  of  this  trace  is  vertically  displaced  only  0.2  m/s  at  the 
scene  from  the  one  used  to  obtain  all  the  rest  of  the  data  but  no  useful  wind 
information  could  be  obtained  from  it. 
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lata  starting  1108  February  14,  1979. 
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The  upper  covariance  function  shown  in  Fig.  17  was  obtained  when  it 
was  snowing  lightly.  Very  likely  the  peak  on  the  function  at  zero  time 
delay  results  from  snow  flakes  going  through  the  field  of  view  near  the 
receiver  since  snow  flakes  farther  away  would  result  in  little  signal  fluctu¬ 
ation  and  even  less  correlation  between  signal  fluctuations  in  the  two 

detectors.  The  remaining  three  covariance  functions  of  Fig.  17  were  obtained 
2  -14 

at  night.  C^  was  4  to  8  x  10  and  apparently  there  was  sufficiently 
different  irradiance  within  the  scene  at  10  ym  wavelength  so  that  performance 
was  similar  to  daytime  operation. 


4.  CONCLUSIONS  AND  RECOMMENDATIONS 

It  is  possible  to  measure  horizontal  crosswinds  by  processing  the  video 
signals  available  in  the  FLIR  thermal  sight.  The  measurement  can  only  be 
made,  however,  under  very  restrictive  conditions.  These  conditions  are  as 
follows: 

1.  C  2  greater  than  about  1  x  10-14  m~2//3 

n 

It  is  unlikely  that  a  passive  wind  measuring  system  can  ever  be  expected 

2 

to  operate  in  all  weather  conditions.  C^  will  vary  at  least  over  four 

orders  of  magnitude  near  the  ground.  While  active  optical  wind  measurement 

2 

systems  almost  never  encounter  so  small  a  value  of  C !  that  operation  is  not 

possible,  a  passive  system  does  because  it  relies  upon  sufficient  ambient 
radiation  from  the  scene,  and  contrast  of  details  within  the  field  of  view. 

In  addition,  passive  wind  sensors  operating  in  the  10  ym  region  have  S/N 
ratios  which  are  about  two  orders  of  magnitude  less  than  a  visual  system. 

2.  Sufficient  contrast  of  the  right  kind  in  the  image 

2 

In  a  passive  wind  sensor,  both  the  reflected  light  and  Cn  tend  to  decrease 
at  the  same  time  and  the  multiplicative  effect  may  result  in  insufficient 
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S/N.  In  the  10  pro  wavelength  region,  a  passive  system  detects  differences 

in  radiation,  which  are  largely  due  to  temperature  differences.  Again, 

2 

Cn  tends  to  be  small  at  the  same  time  that  the  temperature  differences 
are  small,  and  the  multiplicative  effect  results  at  times  in  insufficient 
S/N.  An  additional  problem  exists  in  that  both  detectors  in  the  focal 
plane  must  see  scenes  having  suitable  and  very  similar  spatial  spectra. 

For  example  a  view  of  the  horizon,  bisecting  the  fields  of  view  of  both 
detectors,  is  ideal.  But  if  the  detectors  observe  just  above  or  below  the 
horizon,  or  if  one  detector  sees  the  horizon  and  the  other  one  does  not,  the 
system  may  not  work  at  all.  In  principle,  one  might  have  a  computer  scan  the 
scene  and  find  those  horizontally  spaced,  adjacent  pairs  of  points  which  both 
display  high  signal  fluctuation  levels.  This  is,  of  course,  what  was  done 
manually  to  obtain  suitable  locations  in  the  focal  plane  from  which  we  ob¬ 
tained  the  experimental  results  shown  in  the  previous  section.  However,  not 
all  locations  selected  in  this  way  are  useful,  presumably  because  the  two 
detectors  see  different  scene  spectra. 

3.  Thermal  sight  stationary 

Motion  of  the  receiver  affects  a  focal  plane  system  to  a  far  greater 
extent  than  passive  systems  which  observe  the  intensity  scintillation  pattern 
at  two  locations  at  the  receiver.  These  systems  have  both  much  larger 
angular  fields  of  view  and  nearly  coincident  fields  of  view,  so  that  both 
detectors  observe  the  same  scene.  In  general,  angular  motion  of  the  larger 
field  of  view  produces  less  signal  modulation,  and  of  course  the  changes  are 
the  same  for  both  apertures.  Both  the  slope  and  the  covariances  have  only  a 
second  order  response  to  simultaneous  fluctuations  in  both  apertures,  mainly 
resulting  in  a  reduction  in  S/N.  For  focal  plane  image  measurements,  the 
FLIR  sight  in  fact  had  a  noticeable  vibration  from  its  cooling  unit  as 
originally  set  up  and  it  was  necessary  to  use  wood  blocks  to  damp  it  out. 
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4.  Thermal  sight  image  in  exact  focus 

All  of  the  experimental  results  are  dependent  upon  the  detectors  being 
placed  exactly  in  the  image  plane.  A  quite  different  calibration  and  wind 
weighting  function  will  be  obtained  if  the  detectors  are  inside  or  outside 
of  focus.  In  fact,  if  they  are  outside  of  focus  the  indicated  direction 
of  the  wind  will  be  reversed.  Since  the  FLIR  sensor  is  focused  manually, 
there  is  no  assurance  of  exact  focus. 

5.  Wind  less  than  4  m/s 

The  sweep  rate  of  the  thermal  sight  inherently  limits  the  wind  speed 
that  can  be  observed.  While  the  limit  was  about  2  m/ s  in  the  present 
experiment,  this  was  due  to  the  tape  recorder  bandwidth;  about  4  m/s  should 
be  possible. 

There  is  an  alternate  way  of  making  the  crosswind  measurement.  By  in¬ 
serting  two  new  detectors  inside  the  focus  of  the  present  thermal  sight 
(probably  using  a  beam  splitter)  it  would  be  possible  to  operate  on  the 
principle  of  separate  receiving  apertures  as  used  in  the  earlier  passive 
systems  operating  at  visual  wavelengths.  A  1  kHz  bandwidth  would  be 
adequate  for  30  m/s  winds,  and  better  S/N  would  result  from  the  reduced 
bandwidth  and  continuous  sampling.  This  arrangement  would  eliminate  the 

focus  problem  and  reduce,  though  not  eliminate,  the  problem  of  image 

2 

motion  and  variation  in  the  scene  spectrum.  The  level  requirement 
will  remain  essentially  the  same. 

In  summary,  it  is  clearly  possible  under  laboratory  conditions  to 

use  the  infrared  image  to  obtain  crosswind  information  day  or  night  when 
2 

Cn  is  high  enough  and  the  scene  is  suitable  (for  example  looking  at  the 
horizon  sharply  defined  against  a  cold  sky).  Exact  focusing  requirements, 
vibration  and  movement  problems,  and  other  unfavorable  characteristics 
inherent  in  using  information  gathered  in  the  image  plane  appear  to  make 


i 
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this  technique  unsatisfactory  for  obtaining  routine  crosswind  information. 
Modifying  the  original  thermal  sight  by  adding  several  detectors  inside 
the  focal  plane  eliminates  the  focus  problem  and  reduces  some  of  the  other 
difficulties.  Even  with  a  stationary  platform,  however,  the  modified 
system  could  not  be  expected  to  operate  all  of  the  time  that  the  thermal 
sight  can  produce  images. 
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